Viscosity from elliptic flow: the path to precision 
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Using viscous relativistic hydrodynamics we show that systematic studies of the impact parameter 
dependence of the eccentricity scaled elliptic flow V2/e can distinguish between different models for 
the calculation of the initial source eccentricity e. This removes the largest present uncertainty in 
the extraction of the specific viscosity of the matter created in relativistic heavy-ion collisions from 
precise elliptic flow measurements. 
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Heavy-ion collisions at the Relativistic Heavy Ion Col- 
lider produce strongly interacting matter at extremely 
high energy densities (a quark-gluon plasma), exhibiting 
hydrodynamic behavior. In non-central collisions the col- 
lective flow is anisotropic. The degree of anisotropy, mea- 
sured by the Fourier coefficients v n of the emitted particle 
distributions in the plane transverse to the beam, is sen- 
sitive to the viscosity of the expanding fireball medium 
0, 0, 0, 0]; the largest anisotropies correspond to flu- 
ids with least viscosity. Measurements of the elliptic 
flow coefficient V2 (which, at midrapidity, dominates all 
other v n ) in 200 A GeV Au+Au collisions have been com- 
pared with viscous relativistic hydrodynamic simulations 
of the fireball matter, yielding an upper limit for the spe- 
cific shear viscosity (i.e. the dimensionless ratio between 
shear viscosity 77 and entropy density s) of 2 < 0.5 [^,@. 
This is about a factor 10 smaller than the minimal val- 
ues (typically found near the liquid-gas transition 
of the corresponding ratio measured in all other known 
(real) liquids Q, with the possible exception of strongly 
interacting systems of ultracold fermionic atoms near the 
unitarity limit Q. In this sense, the quark-gluon plasma 
(QGP) appears to be the most perfect liquid ever ob- 
served. 

The uncertainty relation places a lower bound on 
the specific shear viscosity [10(, and explicit compu- 
tation in a large class of strongly coupled field theo- 
ries, using the AdS/CFT correspondence, puts it near 
(s)kss — 17 ~ 0.08 [U]. The above empirical bound is 
sufficiently close to this fundamental limit to have gener- 
ated widespread interest in a precise measurement of the 
specific shear viscosity of the matter created at RHIC. 
Clearly, a measured value that saturates the KSS bound 
would have broad ramifications. It thus came as an un- 
welcome surprise when it was realized that our insuffi- 
ciently precise knowledge of the initial fireball eccentric- 
(which depends on rapidity and, 
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lty e = 

through anisotropic pressure gradients, drives the flow 
anisotropy) introduces a large, apparently irreducible un- 
certainty into the extraction of the specific shear viscosity 
that, for 2 „ O (jh), can be 100% or more @. 

Let us explain the situation in more detail. To date 
elliptic flow appears to be the observable that shows the 



strongest sensitivity to shear viscosity. Two main man- 
ifestations of shear viscous effects have been identified: 
(i) Shear viscosity reduces the amount of elliptic flow 
below the value generated in an ideal fluid [3, 0, 0, 0] ; 
for fixed initial conditions, in particular for a given ini- 
tial fireball eccentricity, the "viscous suppression factor" 
grows monotonically with rj/s [3, 0,01 • (ii) When scaled 
by the initial eccentricity, the elliptic flow for fixed rj/s is 
reduced more strongly in smaller collision systems that 
create shorter-lived fireballs than for larger, longer-lived 
fireballs. Viscous effects are stronger in peripheral than 
in central collisions, larger in Cu+Cu than in Au+Au 
collisions, and they become weaker at higher collision en- 
ergies 

Based on observation (i), Luzum and Romatschke Q 
made a first attempt to extract rj/s from experimental 
elliptic flow measurements in minimum bias Au+Au col- 
lisions at y/s = 200 A GeV. A key ingredient of such an 
analysis is the ideal fluid dynamical baseline with which 
the data are compared to establish the "viscous suppres- 
sion factor" that is used to extract rj/s. In ideal fluid 
dynamics the elliptic flow is directly proportional to the 
initial fireball eccentricity, but the latter cannot be ex- 
perimentally accessed because there are no known probes 
of the reaction zone that escape directly from the fireball 
and probe only the initial state, without any contribu- 
tions from later stages of the expansion. It must therefore 
be estimated theoretically from the overlap geometry cor- 
responding to the impact parameter of the collision which 
can be extracted from measurements of the final hadron 
multiplicity and transverse energy. Unfortunately, the- 
oretical models used to calculate the initial energy and 
entropy density distributions for given impact parameters 
differ by up to 30% in the predicted source eccentricity 
In the study Q two models were studied whose 
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eccentricities e differed by about 20%. The correspond- 
ing 20% variation in the ideal fluid baseline for the elliptic 
flow V2 led to variations by more than a factor 2 in the ex- 
tracted value of T]/s. While the hydrodynamical calcula- 
tions in [B[ made several other assumptions (in particular 
in the equation of state) that affect the ideal fluid base- 
line for V2, their combined effects are likely smaller than 
that produced by the uncertainty in e. More importantly, 
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however, it is known how to systematically improve on 
these approximations in the future and thus dramatically 
reduce their contribution to the systematic error of r]/s. 
On the other hand, it appears impossible to eliminate the 
uncertainty in the source eccentricity on purely theoret- 
ical grounds, and it is not known whether it can ever be 
measured experimentally. With the method used in [Bj], 
a precise meaurement of rj/s, with errors below about a 
factor 2, thus appears to be impossible. 

In this Letter we show that, by exploiting the observa- 
tion (ii) above, one can make measurements that permit 
to clearly distinguish between the two models for e stud- 
ied in [5| . This clears the path to a precise extraction of 
rj/s (or at least of a combination of the specific shear and 
bulk viscosities [(j. fl5|) from elliptic flow measurements. 

Figure Q] shows the initial spatial fireball eccentric- 
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ity calculated from the Glauber and fKLN 
17[ models for 200 A GeV Au+Au collisions as a func- 
tion of impact parameter b. It is defined as e(b) = 

[fcAiZ^X^) where < x ^ b ) is the ener §y den - 

sity in the transverse plane at z — 0; here z denotes the 
beam direction, x the direction of the impact parameter, 
and y points orthogonal to the reaction plane. For the 
Glauber model e(x±; b) is taken to be proportional to a 
superposition of wounded nucleon (85%) and binary col- 
lision (15%) densities [HI]. In the fKLN model the shape 
of e(x±;b) is controlled by the dependence of the gluon 
saturation momentum Q s on the transverse position x±. 
We compute it according to R ef. J 121 . using a program 
kindly provided by the authors [18| . 




FIG. 1: (Color online) Impact parameter dependence of the 
initial fireball eccentricity e (see text for details). 

A quick look at Fig.[T]shows that the eccentricities from 
the fKLN model are O(20— 30%) larger than those from 
the Glauber model (TT|, [T3|. Closer inspection reveals, 
however, that the excess depends strongly on impact pa- 
rameter (see inset): at small b the fKLN eccentricity is 
more than 60% larger than the Glauber one whereas at 



large b > 8 fm the excess drops to < 20%. Had e G turned 
out to be simply proportional to e fKLN , the constant of 
proportionality (and with it the model dependence of the 
initial eccentricity) could have been simply eliminated by 
forming the double ratio (t^/^pcriphcrai/^/^ccntrai and 
exploiting the system size dependence of viscous effects 
to extract rj/s. The inset in Fig. [T] shows that this will 
not work. 

To make progress, let us next look at the sensitivity of 
the eccentricity-scaled elliptic flow v 2 /e to rj/s and then 
proceed to find a way to determine which model for s 
should be used for scaling the experimentally measured 
elliptic flow. The following results are obtained from 
viscous hydrodynamic simulations of 200 A GeV Au+Au 
collisions, with constant rj/s and standard initial and fi- 
nal conditions [1, - We comment on possible effects 
from a temperature dependence of rj/s at the end. 

In Fig. [5] we show w™ b /e m b for minimum bias colli- 
sions, as a function of rj/s. v^ h is obtained from the 
minimum bias pion spectrum dN™ h / (dycPpT) (without 
resonance decay feeddown), and e m b is computed as 
above from the minimum bias energy density e m b(£Cj_) = 



2 f max dbbe(xr,b)/b1 1 ^ 1 with 6 max : 



; 13 fin. 



0.6- 

I ' 5 " 
e^0.4- 

0.3 

0.2 




O Glauber 
OfKLN 



EOSI~**=@._0 0.04 0.08 0.12 0.16 0.2 0.24 
T|/s 




0.04 0.08 0.12 0.16 0.2 0.24 

T|/s 



FIG. 2: (Color online) Scaled elliptic flow vi/e for minimum 
bias 200 A GeV Au+Au collisions, as a function of specific 
entropy rj/s, from viscous hydrodynamics with two different 
equations of state, for the Glauber and fKLN initial state 
models. Inset: The fractional viscous suppression of w"'/^ 
as a function of rj/s. See text for discussion. 

The figure shows that for minimum bias collisions ec- 
centricity scaling works in viscous hydrodynamics, i.e. 
one obtains almost identical curves for different initial 
eccentricity models (the same does not hold at fixed im- 
pact parameters, see below). For any given viscosity rj/s, 
the scaled elliptic flow depends only on the stiffness on 
the equation of state: EOS I, which describes a massless 
parton gas with e = 3p and sound speed c s = l/\/3, 
gives more elliptic flow per eccentricity than the softer 
SM-EOS Q [l9( which matches a massless parton gas to 
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a hadron resonance gas at T c = 164 MeV through a first- 
order phase transition, with c s = in the mixed phase 
and c s w l/v6 in the hadronic phase. However, the frac- 
tional suppression of v^le by shear viscosity below its 
ideal fluid value is almost independent of the EOS: defin- 

(»2° b /E»b) v]: 



ing /„ 



the fraction of scaled elliptic 



flow generated in viscous hydrodynamics relative to the 
ideal fluid value, the inset in Fig. [5] shows that this frac- 
tion is an approximately universal function of the shear 
viscosity rj/s, with practically no sensitivity to the ini- 
tial eccentricity and only weak dependence on the EOS. 
(The somewhat lower fractions for EOS I are caused by 
earlier decoupling in this highly explosive case [4| which 
cuts off the evolution of vi before full saturation [lij].) 
Since the EOS can in principle be obtained from Lattice 
QCD with arbitrary precision, the EOS dependence is 
not a concern. What counts is that there is no large de- 
pendence on the model for the initial eccentricity e which 
can not be reliably calculated from first principles. 

Figure [2] states that if we know w™ b /e m b and the EOS, 
we can determine rj/s. But experimentally one can only 
measure v™ h while e m b must be calculated from a model. 
How can we tell which model for e is correct? The mini- 
mum bias eccentricities for the fKLN and Glauber mod- 
els (e^ LN = 0.197 vs. e^ b = 0.174, see horizontal lines in 
Fig. [T]) differ by 12%, and Fig. shows that the corre- 
sponding ~ 12% uncertainty in v-x/s leads to a factor ~ 2 
uncertainty in rj/s when Tj/s = (9(1/4tt). 



and SM-EOS Q (Fig.©. The left panels in these Figures 
showt;f <LN (6)/e f KLN(&) (top) and vf(b)/e G (b) (bottom), 
respectively. Comparing the upper and lower left pan- 
els in Figs. [3]or[¥]one observes approximate eccentricity 
scaling at all impact parameters, although not with the 
same degree of precision as for the minimum bias average 
shown in Fig. [2] At each impact parameter, the viscous 
suppression of V2/S is grows monotonically with rj/s. 

For EOS I (Fig. [3]), v^js decreases monotonically with 
b, reflecting the earlier freeze-out and a lower degree of 
saturation of the elliptic flow in more peripheral colli- 
sions. For SM-EOS Q we observe a more complex pat- 
tern: as seen in the left panels of Fig. 31 u 2 /£ increases 
with b for the ideal fluid (except for very large b where 
early freeze-out again takes its toll) but decreases with b 
for viscous fluids once rj/s exceeds about once or twice 
the KSS bound (depending on whether we use Glauber 
or fKLN initial conditions). The increase with b seen for 
the ideal fluid is well-known [3] and reflects the effec- 
tive stiffening of the EOS (i.e. a larger effective speed 
of sound) as the system evolves out of the very soft 
mixed phase (which in central Au+Au collisions at top 
RHIC energies suppresses the buildup of elliptic flow) 
into the significantly harder hadronic phase (which dom- 
inates elliptic flow buildup in the more peripheral colli- 
sions). Shear viscosity effectively smoothes the EOS in 
the transition region from a first order phase transition 
into a smooth crossover jij, restoring the monotonic de- 
crease of V2(b)/e(b) seen also in the left panels of Fig. [3] 
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FIG. 3: (Color online) Scaled elliptic flow (v2/e)(b) from 
viscous hydrodynamics with EOS I and either Glauber or 
fKLN initial conditions. The lines correspond to Tj/s = 
0.0, 0.04, 0.24 (from top to bottom in steps of 0.04). 
The left panels show «2 KLN /e£kln (top) and v§ /sg (bottom), 
while the right panels show the "swapped ratios" «2 KLN /eg 
(top) and v§ /eozln (bottom). See text for discussion. 



The solution to this problem is given in Figs. [3] and [4] 
which, instead of the minimum bias ratios, explore the 
impact parameter dependence of v%/e for EOS I (Fig. [3]) 



FIG. 4: (Color online) Same as Fig. [3j but for the equation 
of state SM-EOS Q with a quark-hadron phase transition. 

The key point of this Letter is, however, made by 
the right panels in Figs. [3] and [U In these we show 
the swapped ratios that one would obtain if Nature 
chose fKLN initial conditions but we as physicists scaled 
the corresponding experimentally measured elliptic flow 
^fKLN i ncorr ectly by dividing by the initial source eccen- 
tricity from the Glauber model (top right panels), or vice 
versa (bottom right panels). The qualitatively different 
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shapes of the curves in the upper and lower right panels of 
Figs.[3]and[4]are a direct reflection of the strong centrality 
dependence of the £qkln / £g eccentricity ratio shown in 
the inset of Fig. Q] It causes the swapped ratio W2 KLN /eg 
in the top right panel to drop much more steeply with 
increasing b than either of the correctly scaled ratios, 
overcoming even the stiffening effects at large b from 
SM-EOS Q in the ideal fluid case. More importantly, 
it causes the other swapped ratio vff /efKLN to increase 
with b from central to mid-peripheral (b ~ 6—8 fm) colli- 
sions. This increase holds even for EOS I, over the entire 
range of rj/s explored here, but it is further strengthened 
at RHIC energies when using SM-EOS Q which effec- 
tively stiffens as b increases. With a more realistic EOS 
that exploits the latest lattice QCD data and replaces 
the first order transition by a smooth crossover [201 ] we 
expect a 6-dependence of v^je that interpolates between 
the shapes in Figs. [3] and [4j 

All available experimental data from Au+Au and 
Cu+Cu collisions at RHIC indicate that V2/S falls mono- 
tonically from central to peripheral collisions, irrespec- 
tive of whether one uses ££kln or £g to scale the mea- 
sured elliptic flow V2 



The lower right 
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panels in Figs. [3] and 0] then appear to exclude the pos- 
sibility that the measured vi arises from flow driven by 
Glauber initial conditions. Furthermore, the left panels 
in Fig. 2] exclude the possibility that the fireball medium 
behaves as an inviscid ideal fluid. We conclude that a 
qualitative comparison of existing data on the centrality 
dependence of the eccentricity-scaled elliptic flow sug- 
gests that collective flow in heavy-ion collisions at RHIC 
is driven by fKLN-like initial conditions, and that the 
fireball evolves as a viscous fluid. Extracting the precise 
value of its viscosity requires a quantitative study that 
goes beyond this analysis. What our work provides, how- 
ever, is the basis for a binary decision tree that allows to 
distinguish between the Glauber and fKLN initialization 
models, thereby eliminating the largest prevailing uncer- 
tainty from such an extraction. 

We close with a word of caution: The present analysis 
assumes that the specific shear viscosity 77/s of the fireball 
medium is independent of collision centrality and thus, 
by implication, independent of temperature. Present the- 
oretical knowledge strongly suggests that rj/s increases 
during hadronization and is significantly larger in the late 
hadronic than in the early QGP phase. If the effective 
77/s (averaged over the expansion history) increases dra- 
matically with 6, on account of the larger role played by 
the hadron phase in the evolution of V2 in peripheral col- 
lisions, it may turn the rise of (v2/s){b) in the lower right 
panel of Fig. |4] into a monotonic decrease, similar to the 
one seen in experiment. For this to happen the effective 
77/s would have to increase from, say, 1/4-7T in central col- 
lisions to above 3/47T at b ~ 7fm. While we believe this 
to be unlikely, only an explicit calculation with a realistic 
model for the temperature dependence of rj/s (and that 



also includes bulk viscosity) will allow one to definitively 
rule out this possibility. 
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